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Based on experimental data on the newly synthesized iron-based superconductors and the rele- 
vant band structure calculations, we propose a minimal two-band BCS-type Hamiltonian with the 
interband Hubbard interaction included. We illustrate that this two-band model is able to capture 
the essential features of unconventional superconductivity and spin density wave (SDW) ordering 
in this family of materials. It is found that bound electron-hole pairs can be condensed to reveal 
the SDW ordering for zero and very small doping, while the superconducting ordering emerges at 
small finite doping, whose pairing symmetry is qualitatively analyzed to be of nodal d-wave. The 
derived analytical formulas not only give out a nearly symmetric phase diagram for electron and 
hole doping, but also is likely able to account for existing main experimental results. Moreover, we 
also derive two important relations for a general two-band model and elaborate how to apply them 
to determine the band width ratio and the effective interband coupling strength from experimental 
data. 
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Since the recent discovery of a new iron-based layered 
superconductor intensive efforts have been focused 
on the nature of superconductivity in this materials both 
experimental! 

oretically 13|, LU, [15|, [16|, [17|, [1 
the well-known copper oxide superconductors, this fam- 
ily of materials exhibit higher critical temperatures, 26K 
in LaOo.gFo.iFeAsjp, 41K in CeOi^F^FeAs Q, 43K in 
SmOi_ x F x FeAs [lj and 52K in PrOo.sgFo.nFeAs [TO], 
as well as 25K in hole doped Lai-^Sr^OFeAs |(|. Very 
recently, a number of preliminary analyses have been 
made to unveil the mystery of superconducting na- 
ture, such as the multiband superconductivity, uncon- 
ventional pairing symmetry, electron-doping and hole- 
doping effects, strong magnetic instability of the nor- 
mal state. Experiments from the specific heat mea- 
surements point-contact tunneling spectroscopy [3], 
and infrared reflectance spectroscopy [7[ provided use- 
ful information. For example, according to the point- 
contact tunneling spectroscopy experiment by Shan et 
al. a remarkable zero-bias conductance peak was 
observed at the (110) interface, indicating the possi- 



be too weak to lead such high critical temperatures [20| . 
The possibility of spin triplet superconductivity was also 



suggested [21 



ble presence of nodal superconductivity. Ou et al. [12 1 
performed angle-integrated photoemission spectroscopy 
measurements and their data provided certain support 
for the existence of SDW ordering and an indication 
of unconventional superconductivity. In the theoreti- 
cal aspect, the nature of unconventional superconduc- 
tivity and the pairing mechanism have also been ex- 
plored preliminarily by several groups based on the den- 
sity functional theory (DFT) and dynamic mean filed 
theory (DMFT) 0, Hi, EE EI El- It was pointed out 
that the electron-phonon interaction in this system may 



In this paper, we propose a minimal two-band BCS- 
type Hamiltonian with an effective interband Hubbard 
interaction term included to model the system. The con- 
struction of our model Hamiltonian is based on band 
structure calculation results and intuitive physical pic- 
tures. Taking into account the main features of fermi sur- 
face for the undoped material calculated from the DFT 
and to capture the essential physics of the superconduc- 
tivity and magnetism in the present system, we adopt a 
minimal version of the Fermi surface on a primary two- 
dimensional square lattice in the Fe-Fe plane: one hole 
band around T and one electron band around M points, 
both crossing the Fermi surface in the undoped case. 
Based on rational physical considerations, we introduce 
an effective interband antiferromagnetic interaction and 
elucidate that the effective intraband antiferromagnetic 
coupling could induce the superconducting pairing with a 
d-wave symmetry. Our main findings are: (i) the normal 
state has an SDW order in the undoped case, while upon 
the charge carrier doping the SDW order drops rapidly 
and the superconducting order emerges; in particular, an- 
alytical results for both SDW and superconducting tran- 
sition temperatures are explicitly presented, and their re- 
spective relations to the SDW and superconducting gaps 
are elaborated in details, with the former agreeing well 
with the existing experiment and the latter being quite 
useful for the future experimental verification; (ii) due to 
the two-band (electron and hole) superconducting nature 
of the material, the transition temperature as a function 
of the effective doping density shows a nearly symmet- 
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ric electron-hole doping dependence, accounting for the 
experimental results; (iii) based on a reasonable analy- 
sis, the two-band superconducting state is expected to 
possess a c£-wave pairing symmetry. Moreover, we ad- 
dress how to enhance superconducting T c in this family 
of iron-based materials and how to verify our two band 
model/theory without any fitting parameter. 

We start from a minimal two-band model, which cap- 
tures the essential physics of the multiband unconven- 
tional superconducting state and SDW ordering, 
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where the two band bare dispersions are respectively 
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1 hereafter) based on the 
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band calculations [13|, [hj, L3JZt uM > where ttt-i = 1/ti and 



777-2 = 1/^2 are the effective masses of the hole and elec- 
tron with 7; 1,2 as the effective nearest neighbor hopping 
integrals in the primary square lattice of sites (with each 
site as a unit cell consisting of two Fe atoms), which is 
rotated by the angle of 7r/4 and is enlarged by a fac- 
tor of \[2 x \/2 with respect to the reduced Fe-Fe lat- 
tice, \x is the chemical potential depending on filling, 
and eg' stands for the band offset, where I = 1,2 rep- 
resents respectively the nearly filled valence band (hole 
band) around T = (0, 0) and the nearly empty conduc- 
tion band (electron band) around M = (tt, 7t), as shown 
schematically in Fig. 1. d ka and Cka are the correspond- 
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FIG. 1: Schematic plot of the bare band dispersions of the va- 
lence (hole) and conduction (electron) bands for the undoped 
parent compound (fx = 0). Eq 1 determines the initial carrier 
density before doping. 

ing electron annihilation operators of bands 1 and 2. For 
the hole (electron) band, we have the density of states 
pi_2 = l/(47rfi.2) with the band width Wh, e = 1/Vi,2- 



£q ' is set to give the carrier density at the undoped case, 
i.e. the hole density of band-1 77^ = 2p\e\p , the electron 
density of band-2 vP e — 2p 2 e^\ and the effective total 
electron number per site is (2 + n Q e — 77°). V^, and V^ 2 , 
are the intraband pairing potentials for the two bands. 
V k y denotes the interband pairing interaction. U a s rep- 
resents the effective interband Hubbard interaction term, 
which will be elaborated in the next paragraph. 

band 1 2 
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FIG. 2: Schematic representation of the origin of the antifer- 
romagnetic interaction. 

We now attempt to give an intuitive physical pic- 
ture to understand the origin of superconducting pair- 
ing and the SDW ordering in this system, ft is well 
known that the superexchange antiferromagnetic cou- 
pling derived from the strong coupling of the single- 
band Hubbard model provides a driving force for the 
unconventional <i-wave superconducting pairing state in 
cuprate superconductors. In the present case with the 
two bands, the 7-th site Hamiltonian is given by Hi = 
Unu a n iU j + U'n ila n atJ , - J H a a ■ cr i2 with rii and cr.; 
as the 7-site electron number and spin operators while U, 
U', and Jh being respectively as the on-site intraband 
Hubbard repulsion, interband one, and the Hund's cou- 
pling constant, which involves also the Hund's coupling 
and interband U' and thus introduces more complica- 
tions. To simply this complicated and subtle issue but 

without loss of the key point, let us consider only the 

(2) 

energy range above the bottom of band 2 (i.e., /i — ), 
below which only the band 1 is available for the occu- 
pation by electrons and thus the Hund's coupling plays 
no role. The number of electrons in this range at zero 
temperature is (1 + ph/ p e )n®N with N as the number of 
sites. For the two empty neighboring sites, each having 
two levels, there will be four levels available for the occu- 
pation by four electrons in the considered energy range. 
In the absence of hopping, due to the on-site Hubbard re- 
pulsion, each site would have two electrons to lower the 
energy (see the left panel of Fig. 2). In this case, there 
are six states, where the two lowest energy states have 
two electrons with the parallel spins in the two different 
bands (both up or down) due the Hund's coupling, while 
the two relatively higher energy degenerate states have 
two opposite spins in the two different bands, with the 
energy difference Jh/2 from the lowest energy level; the 
rest two states have an even higher energy with opposite 
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spins in either of two bands as the intraband on-site U - 
repulsion is stronger than the interband one ([/') (Here 
we refer to these six states as being in the lower Hub- 
bard band). When the hopping is included, more very 
high energy states, with each site having at least three 
electrons, are involved (we refer to these states as be- 
ing in the upper Hubbard band). If an electron hops 
from the bare ground state via virtual processes to the 
upper Hubbard band and then back, which involves an 
additional intraband U and corresponds actually to the 
antiferromagnetic superexchange-like process, lowers the 
total ground state energy of the two sites. A typical hop- 
ping process is schematically shown in Fig. 1. Analogous 
to the single-band case, the present superexchange-like 
interaction could also lead to the effective antiferromag- 
netic interaction between the two neighboring sites and 
thus antiferromanetic fluctuations in the whole lattice (as 
only a part of mobile electrons/sites are involved), which 
introduce the effective superconducting pairing interac- 
tion having most likely a nodal d-wave symmetry, as in 
cuprates. Also for simplicity but with out loss of main 
physics, here we wish to argue that only the interband 
term is retained [23| in the effective Hamiltonian (1) af- 
ter the virtual processes because the intraband U may be 
significantly greater than the interband V. [HI Since the 
interband repulsion term lifts the energy Jh/2 above the 
bare ground state energy level, U e g is ~ Jh/2. Below we 
elaborate first that this C/ e ff-term can indeed lead to an 
SDW order at zero or very small electron/hole doping, 
where the superconducting paring order is suppressed. 

The SDW transition temperature Tsdw can be calcu- 
lated from the following equation [23| 



i = £/cffx£ 2 (Q), 

where Xq 2 is the interband spin susceptibility 



S 2 (Q) = E 



/(£ifc) ~ /C&fc+g) 

£,lk — £,2k+Q 



(2) 



(3) 



with /(e) as the Fermi-Dirac function f(e) = 1/(1 + 
e £ / Tc ). To obtain a simple analytic formula of TsdWj 
here we set m\ — vn-i and eg = = eoj where the 
perfect nesting with vector Q = (tt, 7r) between the two 
bands occurs at the undoped case(/z = 0), which leads to 
the SDW instability. Integrating the RHS of Eq. we 
have an equation for Tsdw 



csdw 
W 



TT V W \ WJ 



8T SDW * 



(4) 

where 7 ss 0.577 is the Euler constant, x is the effective 
small doping, and the condition eo 3> Tsdw has been 
used. When the electron or hole doping is zero, i.e. x = 
0, we have the largest Tsdw- Remarkably, Tsdw drops 
drastically whenever more electrons or holes are doped, 



as seen clearly in Eq.Q and Fig. 3. If mi 7^ mi, Tsdw 
is expected to be lowered. 

Below Tsdw, the SDW ordering emerges, whose order 
parameter may be defined as 

t/eff 



A 



SDW 



( CfeCT 4 



k+QS 



(5) 



where (• • • ) denotes thermodynamic average. Asdw sat- 
isfies the following equation 

1 tt ffok + tok) - fimk- ^) ,„v 

1 = t/eff 2^ ^ ' ( 6 ) 



where flk ■ 

V2k = (Clfc 



Vvlk + A sdw Vik = (Cifc - 6fc+Q)/2 and 
^2fc+Q)/2. Here the effective repulsive in- 
teraction (C/eff ) between interband electrons may also be 
viewed as an attractive pairing interaction between elec- 
tron and hole. Physically, it is noted that the pairing 
leads to form a " condensate" of bound electron- hole pairs 
in the triplet state or "excitons" (SB)], which exhibits the 
SDW ordering. The condensate of electron-hole pairs is 
actually an counterpart of Cooper electron-electron pairs. 
In this sense, it is straightforward to obtain the famous 
relation 2Asdw(0) ~ 3.5Tsdw, as in the case of the con- 
ventional weak-coupling superconductivity. Here we pin- 
point out that this relation can be quantitatively verified 
by independent experimental methods, e.g., the optical 
conductivity spectra Q and resistivity or specific heat 
measurements. As seen from Ref.[8|, it may be estimated 
that 2A SD w(8K) « 350cm- 1 = 504K and T SD w ~ 150K, 
leading to 2Asdw(8K)/Tsdw ~ 3.4, which agrees well 
with the present theory. It is worthwhile noting that, 
since the effective hopping integrals tip in the present 
model correspond to those between nearest neighbor- 
ing sites (unit cells) in the primary lattice, the above 
SDW ordering pattern is stripe-like antiferromagnetic on 
the reduced Fe-Fe square lattice, namely, the spin or- 
dering pattern is ferromagnetic in each stripe along the 
x'(or y')-direction of the reduced Fe-Fe lattice, while 
it is antiferromagnetic between stripes in the y'(or x')- 
direction. This kind of SDW ordering was likely observed 
in a very recent neutron scattering experiment (26j . Re- 
markably, we estimate from Eq. ([5]) that the antiferro- 
magnetic moment per Fe-atom (in unit of Bohr mag- 
neton /i B ): m A F ~ (4/2) x m~ = 2 J2kcr( c ka d l +Q9 ) = 
4A SD w/E4ff = 2 x 3.52 (T c /W)/(U cS /W) sa 0.31, where 
mi = {—iy J2a( c ivd]g) is the i-th site moment [27]. This 
estimation is also in agreement with the data reported in 
Ref [26j . In addition, we may have one more expecta- 
tion from Eq.([J| that Tsdw is normally decreased with 
the increase of pressure because the effective band width 
W (or the effective hopping integral t) that dominates 
Tsdw m the exponential term is increased, as seen in 
Ref. [H. 

At this stage, we turn to address the superconduct- 
ing ordering. In view of the fact that the SDW or- 
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der drops very sharply to zero at a very small critical 
doping, it is reasonable and convenient to ignore the 
effect of SDW order on the superconducting ordering 
above the critical doping level. In the following calcu- 
lations, the pairing potentials involving two bands are 
expressed as: Vlh = Jhhlklk* , V3 = Jeelklk' and 



V kk' 21 = Jhe,eh~fkjk>, where J hh , J ee , and J he ,eh are the 



corresponding coupling constants, jk — cos(L8k) with 
L = 2, 1, denoting respectively the nodal d-wave, nodal 
p-wave, and isotropic s-wave pairing functions, and Ok 
being the angle between the vector k and the x-axis. 
Taking the BCS mean field approximation, the quasi- 
particle eigenspectrum of the Z-th band is then given by 
Eik = y/£,f k + |A;| 2 7^, where |A/| is the gap amplitude 
of the l-th band(/ = h, e). They are determined from the 
following coupled gap equations 29]: 

A/, = ^7fe(Jftft(c_ fc |C feT ) + J he (d-kidk^}), 

k 

A e = ^7fc(J ee (d_ fci d fcT ) + J e ft(c_fc^c fcT )). (7) 
fc 

Then the self-consistent gap equations for Aft and A e 
read 



J hh K x J he K 2 W Aft 
Jeftifi J ee K 2 \ A e 



A/, 
A P 



(8) 



where K h2 = £ fc 7^ tanh (£ li2fe /2T) /E h2k satisfy the 
following equation 
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The superconducting transition temperature T c is ac- 
tually determined from Eq. ([9]) with A; — > 0, which de- 
pends on the pairing interaction strengths and the normal 
state dispersions of the two bands. 

Let us first address a more general case: J ee , Jhh > 
0, and J ee Jhh ^ J eh J he > (noting that J e h = Ji e )- 
By the introduction of dimensionless couplings: Jhh,ee = 

Jhh,ee/Wh,e, JehMe = Jeh,he/Wh,e, and J J = J e hJhe 

J, , Ji . . the T c -formula is derived as 



VWeWft 



2 Cr <= [n e n h {2 - n e ){2 - n h )Y e 



(10) 



with C Tc = ln( e T/7r) w -0.568 (for L = 2, 1, cases) and 
A ro d as the reduced pairing strength being given by 



red 



1~ n e (2-ne)W? Jhh-Je. 

— J J ill ; ^ n H 

4 nft(2 - n h )Wl 2 



Jeh J ht 



+ Jhh) > /JJ 



where n e = n° e + [W h /(W e + W h )]x and n h = n° h - 
[We/(W e + Wh)]x, and the condition e{, 1,2) > T c has been 



FIG. 3: (Color online) Superconducting critical temperature 
T c and the SDW transition temperature Tsdw as functions 
of the effective doping x. The Tsdw curve is calculated from 
Eq. flU with the parameters U e g/W = 0.3 and e /W = 0.05; 
the T c curves are plotted by simply substituting the relevant 
reduced parameters into Eq. (|10|l . where J ee /W = 0.15, and 
| Jeh\/Jee = 0.3, 0.6, 0.8, respectively for Jhh = Jee and Jhh = 
0.8J ee . 



used. For a special case: J J = 0, by taking the limit of 
J J —>■ in the above equation, the T c -formula is reduced 
to 



VWeWh 77 




y/n e {2 - n e ) 



(11) 



In this special case, more intriguingly, it is found from 
Eq. ([5]) that the ratio between two gaps r(T) = Aft/A e — 
Jhe/Jee = Jhh/ Jeh, being independent of temperature 
and other variables. Therefore, if the ratio value \r\ is 
experimentally found to be independent of temperature 
and other variables, it is implied that the system is just 
in this special case. Note that this result is valid for any 
two-band superconductivity model described by Eq. ([7]). 
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For a general case, i.e., J J ^ 0, the phase diagram of 
T c -x calculated from above formula is plotted in Fig. 3, 
where we choose the parameters as W e — Wh — W (for 
simplicity but without loss of generality) , J ee /W = 0.15, 
and e^ 1,2) = e = 0.05W, U cS /W = 0.3 (since J H « 0.9 
eV, i.e., U c ff ~ 0.45 eV, such choice leads to the effec- 
tive band width W ~ 1.5 eV, which is not unreason- 
able in the present system). From Fig. 3, we note that 
the superconducting order emerges from a very small ef- 
fective doping and then decays rather slowly to zero at 



:(e,h) 



When J ee = Jhh and n h = 



the T c vs. x dependence is symmetric for the effective 
electron and hole doping ; while the symmetric feature 
changes slightly if J ee (or n°) is different from Jhh (or nP e 
) not too much (as estimated from the band calculations), 
in agreement with the experimental result Q . In partic- 
ular, if one intraband pairing strength (or the interband 
coupling) with the largest value is fixed, both the inter- 
band coupling (regardless of its sign) and other intraband 
pairing strength(s) enhance T c significantly, reaching the 
maximum as \J e h\ ~ > Jhh ~* Jee, whose value is explic- 
itly given by Eq. (fTT|) . This feature may be helpful for 
searching even higher T c superconductors. Based on our 
results, another possible way to increase T c in this fam- 
ily of superconductors is to increase the effective band 
width W (or the effective hopping integral t ) since we 
have approximately T c cx We~ a ^ w with a being a W- 
independent coefficient by noting that J e e,hh/W cx W. In 
this sense, a higher pressure may enhance T c , in contrast 
to that of Tsdw- This expectation was seen in a very 
recent high pressure experiment that showed clearly the 
enhancement of T c with the shrinkage of the lattice [3(| • 
The normalized tunneling conductances along the di- 
rections (110) and (001) of the primary lattice versus the 
bias voltage are calculated self-consistently from the gap 
equations ([8]) with the nodal d-wave (cos(20fc)) pairing 
at zero temperature and electron doping at x — 0.1, as 
shown in Fig. 4 for a set of parameters used in plotting 
Fig 3 (except for the upper panel where Jhh — 0.5J ee ). 
A sharp zero-bias-conductance-peak (ZBCP) along (110) 
and two coherence peaks corresponding to the two gaps 
(a larger A° and a smaller A?) are clearly seen, as ex- 
pected. In addition, the two weak kinks in the ZBCP 
curve can also be seen, which correspond to the two gaps 
as well. The present ZBCP results are in agreement with 
the experimental observation reported in Ref. Q 

Finally, we wish to address how to verify unambigu- 
ously the present two band model/theory experimentally, 
without any fitting parameter. From the definitions of 
K\,2 below Eq. ([7j), we have 



K li2 (T c 
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FIG. 4: (Color online) Normalized tunneling conductance as 
a function of bias voltage at zero temperature for an electron 
doping at x — 0.1. The parameters are Wh = W e = 1.5 eV, 
Jee/W = 0.15, \Jeh\/Jee = 0.3 and J hh / Jee = 0.8 for the 
lower panel, which gives rise to the self-consistent A° = 3.93 
meV and — 2.75 meV; for the upper panel Jhh/ Jee = 0.5, 
which leads to A° = 3.06 meV and A° = 1.46 meV. 



where Co = (In 2 — 1/2) sa 0.193 for the nodal d-wave 
and p-wave cases (while it is zero for the gapped s-wave 
and p-wave cases). On the other hand, from Eq. (|8|), we 
have W h ,eK ia {T) = [-J ee , hh + r- l > l {T)J he>eh }/Tj for 
any temperature T < T c . Combining them, we can find 
an important relation: 

ln(|A°|/r c ) + Co-C T J . ,W h 

x \r Q r c \ - —. 



Co — Ct c 



(12) 



|-ln(|A°|/T c ) 

where ro iC = r(0),r(T c ) are respectively the above intro- 
duced gap ratio r(T) below Eq. (JTTJ) at zero and transi- 
tion temperatures. Since the RHS of Eq. (fi"2l) depends 
only on the ratio of the two band widths, this relation can 
be verified by experimentally measured data for |Ag fe | 
and \ro tC \ at various doping levels and then be used to 
determine the ratio of Wh/W e . Similarly, we can obtain 
another useful relation: 



ln|r | 



\Jeh\ 



(13) 



IT / 



h . e 



[l + (W h /W e )\r r c \-i](\r c \-\r \) jj ' 

Once (Wh/W e ) is determined from Eq. (fT2")) . this rela- 
tion can not only be checked by experimental data of 
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ro, c but also be used to determine the effective interband 
coupling strength \ J eh \/JJ. It is also interesting to see 
from Eq. (fP5j) that (i) as long as |r | w \r c \ is observed 
experimentally, the system is approximately in the men- 
tioned special case; (ii) if |ro| ~ 1 is seen experimentally, 
except the case (i), the interband coupling is negligible; 
(iii) if two bands superconduct independently with two 
different T c 's, the system T c is determined by the band 
having a higher one, and the LHS of Eq. (fl3|) approaches 
to zero because r c is either divergent or zero at the lower 
T c . It is notable that the above results are also valid for 
the two-band s-wave and gapped p-wave superconduc- 
tors similar to the present system, except for a different 
constant term (Co — Ct c ) s , p ~ 0.568 . 

In summary, we have proposed for the first time a min- 
imal two-band (hole and electron) model, with the inter- 
band Hubbard interaction being also included. It has 
been shown that this simplified model is able to capture 
the essential physics of unconventional superconductiv- 
ity and spin density wave ordering in the addressed new 
family of materials. The present theory not only accounts 
for very recent experimental results, including the emer- 
gence of hole-doped superconducting states and the neu- 
tron scattering analysis on the SDW ordering, but also 
demonstrates/elaborates a key role of the interband pair- 
ing coupling played in the significant enhancement of T c 
(in comparison with that of each band) for a generic two- 
band system, regardless of the pairing symmetry. 
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